A HYPER-GEOMETRIC APPROACH TO THE 
BMV-CONJECTURE 
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Abstract. We provide a representation of the (signed) BMV-measure by sto- 
chastic means and prove positivity of the respective measures in dimension 
d = 3 in several non-trivial cases by combinatorial methods. 

1. Introduction and Results 

We aim to provide a representation of the (signed) measure related to the Bessis- 
Moussa-Villani conjecture (in the sequel BMV) by stochastic methods and calculate 
non-trivial cases in dimension 3 by hyper-geometric methods. 

Definition 1. Let d > 1 be fixed. Let A, B he complex, hermitian d x d matrices 
and B > 0, then we denote 

0^'-^(z) :=tr(exp(A-zB)) 

for zeC. 

The Bessis-Moussa-Villani conjecture (open since 1975, see j^Jj asserts that the 
function cj)^'^ is completely monotone, i.e. cj)"^'^ is the Laplace transform of a 
positive measure fi^'^ supported by [0,oo[, 

tr(exp(A - zB)) = / exp{~zx)n^'^ (dx). 
Jo 

Since the function (f)^'^ is always Laplace transform of a possibly signed measure 
on [0,00 [, we shall always denote this signed measure by fj,^'^ . 

The BMV conjecture is closely related to convergence assertions on perturbation 
series in quantum mechanics and there is a substantial literature on it (recently jSj 
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has been published, where several further references can be found). We quote from 
"The BMV conjecture would entail a number of interesting inequalities not 
just for quantum partition functions, but also for their derivatives (a badly needed 
tool). Despite a lot of work, some by prominent mathematical physicists, only 
some simple cases have been decided. So far all results, including fairly extensive 
numerical experiments, are in agreement with the conjecture" . As an example 
of recent progress on positive results we mention that the BMV-conjecture was 
shown to hold true in an average sense in ^l- Originally the BMV-conjecture 
was formulated more generally, namely, that z i— > (e,exp(A — zB)e) is completely 
monotone for each eigenvector e of B. This first conjecture was seen to be wrong 
immediately (see the end of the article 0). In Appendix 2 we provide a simple 
counter-example for the sake of completeness. 
At first sight the following relations hold true: 

Proposition 1. Let A,B be hermitian dx d matrices, B > 0, then (j)^'^{z) > 
for z > and 

^4>^'^ (z) = - tr(exp(yl - zB)B) 
A-^-^.Bj-^-) ^ exp(-s(A- zB))Bexp{s(A - zB))Bds exp{A - zB)^ 

for z>0. Hence -£(t)^'^{z) > and -£^(j)^^^{z) >Oforz>0. 

Proof. The first assertion follows from the fact that the eigenvalues of exp(^ — zB) 
are non-neagtive and the second from the derivative of the function exp off (see 
for instance 7 , Theorem 38.2), 

d 

— exp{A - zB) = - exp{A ~ zB) I exp(-s(yl - zB))B exp(s(A - zB))ds. 
dz Jq 



Hence 

d 
dz 



(j)^'^{z) = -tr (^cxp{A- zB) cxp{-s{A- zB))Bexp{s{A- zB))ds 

= - [ tr {exp{A - zB) exp(-s(A - zB))B exp(s(A ~ zB))) ds 
Jo 

= - / tr(exp(A - zB)B)ds 
Jo 



= -tT{cxp{A- zB)B). 
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The second formula follows by a similar reasoning. We conclude the inequalities in 
a "moving frame" associated to the eigenbasis of exp s{A — zB). □ 

Bernstein's Theorem (see for instance jS]) tells that a smooth function (f) : R>o 
M is the Laplace-transform of a non-negative measure /i on M>o if and only if 

1)"(/)*^"^(2:) > for z > 0. For the BMV-function (j)^'^ we know by the previous 
Lemma at least, that the Bernstein condition holds for n = 0, 1,2. For dimension 

> 3 the case ti = 3 is unknown in general. Having Bernstein's Theorem in 
mind, we see that the validity of the BMV-conjecture is equivalent to a sequence 
of interesting trace inequalities for hermitian matrices. 

The following simple transformation properties are immediately proved. 

— T — T 

(1) Given a unitary matrix U in dimension d, then yP^^ ,ubu _ ^a,b ^ This 
is due to the unitary invariance of the trace functional. 

(2) Let Id denote the identity matrix in dimension d. Then ^-^-fAi/^.B _ 
exp(Ai)/i"*'^ for all real Ai, since the identity matrix commutes with A, B. 

(3) ^j_A^B+X2id — -f A2) for A2 > — 6min, where 6min denotes the minimal 
eigenvalue of -B, since a translation of B by X2ld corresponds to a translation 
of the measure by A2. 

Furthermore the following cases are known, where the BMV-conjecture holds 
true. 

(1) If A and B commute, the BMV-conjecture holds true. 

(2) If d = 1, 2, the BMV-conjecture holds true. 

(3) If B has at most two different eigenvalues, the BMV-conjecture holds true. 

(4) Let _B be a diagonal matrix. If the off-diagonal elements of A are non- 
negative, the Dyson expansion (see Section 2 for a stochastic proof) yields 
that the BMV-conjecture holds true. 

In view of all these well-known facts (for more investigations in these directions 
see ^), the first non-trivial case, which appears in lowest non-trivial dimension, is 
the following. Take d = 3, B = diag{bi,b2,b^) a diagonal matrix and A = (aij), 
and assume 012013023 < 0. In this article we give - by hyper-geometric methods - 
a partial positive answer in this case. 



4 MICHAEL DRMOTA, WALTER SCHACHERMAYER, JOSEF TEICHMANN 

We first provide a Feynman-Kac-type construction - by stochastic means ~ of 
the Dyson series. Since Feynman-Kac-type Theorems are often used to prove that 
a function is a Laplace transform, we were motivated to construct an appropriate 
Markov process for the non-stochastic matrix A in order to prove a Feynman- 
Kac representation for (p^'^ . From this representation we can deduce the Dyson 
expansion and we are able to deduce a "semantics" of the problem, however, we were 
not able to conclude the result directly. The proof of the Feynman-Kac Theorem can 
be found in the Appendix 1, its application in order to prove the Dyson expansion 
in Section 2, see Theorem El 

In Section 3 we concentrate on the 3-dimensional case, where we meet an impor- 
tant combinatorial simplification, see 12. 71 which then leads to a summation problem 
in the theorie of hyper-geometric series. Finally we are able to prove the following 
result. 

Theorem 1. Given a real, symmetric 3x3 matrix A = (aij) and a diagonal matrix 
B = diag(hi,b2,b^) with diagonal elements < 6i < 63 < 62. We assume that the 
following two conditions hold true: 

(1) kial |qi2| |a23| 

(2) 011(62 - h) + a22(&3 - bi) + 033(61 - 62) > 0. 

Then the function (j)'^'^{z) := tr(exp(j4 — zBj) is completely monotone and the 
BMV- conjecture holds. Furthermore the BMV-conjecture holds (trivially) if two of 
the three eigenvalues 61,62,63 agree or 012013023 > 0. 

Remark 1. The unusual order 61 < 63 < 62 is due to the structure of our proof, 
see Section 3. Later we shall assume 61 = 0, which is possible without restriction 
of generality as we have noted above under transformation property (3) on page 3. 

Remark 2. The two conditions in (1) are related to positivity on the intervals ]0, 63[ 
and ]63,62[, respectively (in this order). The second condition is a linear functional 
on the diagonal values of A and appears to be the same on both intervals. 

Remark 3. The proof of Theorem^will be given in Section 5. We note, however, 
that the assertions of the last sentence can be proved immediately (the argument 
shows all trivial cases). 
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Proof. Assume that 012013023 > 0, then we can make a change of coordinates such 
that aij > 0, for « 7^ J by multiplying two coordinates by —1. If aij > holds for 
i ^ J, then by Theorem |21 the measure fi^'^ is a sum of non-negative measures, 
hence non-negative. If 62 = ^'3, then B has a 2-dimensional eigenspace, where we 
can rotate without changing B, consequently we can find an orthogonal matrix U 
such that {U'^AU)23 = and U^BU = B. The trace is invariant under rotations, 
so 

(jj'^'^iz) = tr(exp(A - zB)) = tr(exp(;7^^C/ - zU'^BU)), 
hence we find ourselves in the first trivial case. □ 

2. Representation of /i^'^ 

In this section we fix d > 2 and a dx d hermitian matrix A. We shall construct a 
Markov process (i^/^'*')o<t<i := iZt'''^\ Xl)o<t<i for (C,«) € S, which leads via the 
Feynman-Kac formula (see Theorem|31in the Appendix 1), to a representation of the 
BMV-measure fi^'^ . The state space of this Markov process is 5* := C x {1, . . . , d}. 
We shall assume a filtered probability space {^t)o<t<i, P), which allows for 

the subsequent constructions. Later we shall identify this space with the polish 
space of cadlag paths on [0, 1] with values in S. 

Let {Nt)t>o be a standard Poisson process with A^o — and jump intensity 
d — 1 defined on and adapted to {^t)o<t<i, P)- We denote by (T,„)„i>o the 

jumping times of (A''()o<t<i, i.e. T,„ := inf{0 < t < 1, Nt > m}, where the infimum 
over the empty set equals infinity. We shall denote by {Xt)t>o the cadlag-process 
starting in a uniformly distributed way at points in S, i.e. P{Xo — i) — ^, and 
having the properties, 

(2.1) T„ = mi{t > r„,_i with Xt ^ Xt^_, } 

(2.2) P(At„ = A;|At„- = I) = 

a — 1 

for m > 1 and fc 7^ / G {1, . . . , rf}, i.e. the process jumps at Poissonian jumping 
times Tm in a uniformly distributed way to another state. Note that this process 
is stationary, i.e. P{Xt = i) = 2 G^idi < i < 1 and 1 < i < d. 
We define for < t < 1 
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where the product is almost surely well defined as A^i < oo almost surely. The 
empty product is defined to be 1. We set 

Yt:= {Zt,Xt) 

for < t < 1. Then (l^)o<t<i := {Zt,Xt)o<t<i is a process with cadlag paths 
starting in a uniformly distributed way at {(1, 1), . . . , (1, d)} in S. 

Remeirk 4. We may and will choose the polish space of cddldg paths on [0, 1] with 
values in S as probability space O with the Borel probability measure P, such that 
the coordinate process on fl together with the canonical filtrations (^t)o<t<i satisfy 
the above requirements. Hence the process is a well defined map on the entire 
probability space O, which will allow us to leave out the usual "almost surely" at 
several occassions. 

We define probability measures P' on {fl,J^, {J-'t)o<t<i) by conditioning on the 
event Xq = i for i e {!,..., rf}, i.e. P' := P{.\Xo = i). With respect to 
the probability measures P' we define, for {C,i) G S, a process (5^^^''')o<t<i := 

o<t<i, P') through 

Xl = Xt 

for < t < 1. 

Proposition 2. Let A be a hermitian d x d matrix. The family of processes 
(y^^^''^)o<t<i on {Ct,J^, (^t)o<t<i, -P') for {C,i) € S defines a Markov process with 
generator 

d 

Af{C,i) = Y,{f{Caij,j)-f{C,i)). 

for all f e C{S,C) and {C,i) G S, where C{S,C) denotes the set of continuous 
functions on S. 
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Proof. Fix / e C{S, C) and (C, i) G S, then 

ii?p.(/(Zp\x,^)-/(C,z)) 

^ oo 

= -^T.^P' (/(^l'^'^ ) - /(C, i)\Nt = k)P{N, = k) 



E(/(c«»^'JV/(c,*)), 



as i ^ 0. □ 

The Feynman-Kac formula allows for a stochastic interpretation of the BMV- 
measure iJ^'^ . Fix r e C^, then Theorem O asserts for functions f^{C, i) '■— riC, for 
(C, i) € S, the following formula to calculate exp{t{A — zB))r. 

Corollary 1. Given r £ C^, a hermitian d x d matrix A and a diagonal matrix B 
with non-negative entries 61, ... , hn, we have 



El^expiJ a{Xs)ds^zj h{X,)ds)f {Yi)\Xn = i 

= Ep. {^mJ^ a{Xl)ds ~z j'^ b{Xl)ds)r{Yl'''^)^ 

= e-^'^-^\exp{A- zB)r)i 

holds true for all z £ C and 1 < i < d. Here a — (ai, . . . , a„) denotes the (real) 
vector of diagonal elements of A. 

This formula allows for an interpretation of z i-^- (exp(yl — zB)r)i as Laplace 
transform of the random variable b{Xs)ds under the (signed) measure on fl 
dO^ 1 

■exp(/ a{Xs)ds)f-{Y^)^x,=^}. 



dP P{Xo ^ i) 
since b{Xs)ds appears linearly in —z in the exponent. 

Given a hermitian d x d matrix A and a diagonal matrix B with non-negative 
diagonal entries bi, . . . ,bd, we define an Jx)-measurable random variable / on C x 
{1, . . . , d} X ($7, JFq, P) to obtain a closed formula for the trace z ^ tr(exp(A — z_B)), 
namely 

d-i^ / 1 if ^0 = « 



(2.3) /(C,z):-de''-^C 
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for C £ C and i = 1, . . . ,d. By Corollary ^ and Definition (j2.3|l we obtain, using 
the notation of Theorem ^ 

d 



i=l 
d 



Y^EMcmI a{Xl)ds~z f 
i=i Jo Jo 

d „i .1 

E{exp{ / a{Xs)ds - z / b{X, 
Jo Jo 

(2.4) =£;(^exp(^ a{Xs)ds-zJ^ b{Xs)ds)f{Y, 



biXl)ds)nY,))e^-'d^ 



.)ds)f{Y,)\Xo=i)P{Xo = t) 



'1 



for z eC. 

We now derive a series representation of the measure fj.^'^ . The function /(Yi) 
can take non-zero values only at loops, i.e. Xq = Xi. First we introduce the subset 
ri„ C consisting of those paths which form loops in {l,...,rf} on [0,1] with 
precisely n jumps for n > 2, i.e. f7„ :— {Xq — Xi,Ni = n}. We define the set 
C„ C {1, . . . , d}" as image set of the path random variable 

p„:f]„^{l,...,dr 

LO (Xti-(cj), . . . ,Xt„-(w)) 

So the subset C„ of {1, . . . , d}" is characterized as set of all n-tuples such that 
no neighbors are equal and the last element Xt„~- is different from the first one 
Xti- = Xq. We denote C := U„>oC„. Elements 7 G C„ are called favorable 
paths of length n. 

The map ord associates to 7 € C„ a monomial in the variables a^ , which is 
called order of the path. The quantities ^^(7) are the respective powers of Uij 
in the monomial ord(7): for 7 e C„ we define 



(2.5) ord(7) a-^-^j^a^^j^ . . ■ o^y^i-iin^inii 

(2.6) =n4 



i<j 



The characteristic char(7) = (fci(7), . . . ,kd(j)) of a path 7 g C„ is defined by 
the number kj (7) of visits in state j 

kj{'-f) := #{Z such that 7; = j}. 
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Notice that the following formula holds for 7 G C„, 
(2.7) i^Z,^.(7) = fc,(7), 

which leads in dimension 2 and 3 to one-to-one relations between char(7) and ord(7) 
(see Lemma 

We shall denote by A„ the n-simplex in M"+^, i.e. the set of vectors (ii, . . . ,in-i-i) S 
with X]r=!i^ ti = 1 and ti > 0. On the n-simplex we shall consider the normal- 
ized uniform law A„, i.e. A„(A„) = 1. For 76 C„, we consider the setp~^(7) C r2„. 
On Pn^{^) we consider the conditional probability := P(-\pn = 7) and the ran- 
dom variable 

dur : p~^(7) A„ 

via dur :~ (Ti, T2 — Ti, . . . , T„ — Tn-i, 1 — T„), which has a uniform distribution A„ 
on the simplex A„ under Pj. Indeed since the conditional distribution of (Ti, T2 — 
Ti,T^ ~ T2, ■ ■ ■ ,Tn — Tn-i) for n > \ under the condition Ni — n is uniform, we 
can conclude the result, i.e. 

P^{Ti e [h,h + dti],T2 -Tie [t2, t2 + dt2], . . . , r„ - T„_i e [i„, i„ + dt„]) 



d n! 

dii...dt„e-('^-i)(d-l)"i^ 



1 e-C'i-i) 



nldti . . . dtn 



for X]r=i — 1- Here we apply that set p„^(7) for 7 G C„ has probability 

since the probability for a trajectory to have n jumps is 

On the simplex A„ we define for a vector h E M"+^ the real- valued random 
variable y'^ — ^i^ij and for g € R"+^ the measure with 

dQs 

i—l 

with respect to the uniform distribution on A„. The image measure of under 
is denoted by r/^"'*, which is a measure on R with support in the convex hull of 
hi, . . . , hn+i- 

Theorem 2. Let A be a hermitian matrix, B a diagonal matrix with non-negative, 
mutually different entries bi,...,bd- The diagonal elements of A are denoted by 
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ai, . . . jflrf. Then the measure (see Definition^ fi^'^ is a signed measure decom- 
posing into an absolutely continuous and singular part 

d 

(2.9) ^i^''^{dx) = exp(a,)(5h. (dx) + ij^'^(x)dx. 

1=1 

ip^'^ is a piecewise continuous function with possible disconituities at bi and with 
support in [minj 6i, max5i]. We have 

(2.10) V^'-^(a;) = ^ (t){j,x)aj-^j^ . . .a^„7i, 



where the density 4'{j, x) is defined by 

1 

n 



(2.11) <j>{j,x)dx := -i,y'^^i>-'''^"'''^i;^i'->^...^i(da;) 



for 7 e C„. 

Proof. We can decompose the Feynman-Kac formula (|2.4|1 by the number of jumps 
iVi, which appear up to time 1. This leads to 

E{exp{f a{X,)ds^z[ b{Xs)ds)f{Yi)) 
Jo Jo 



oo „i „i 

E{exp{ / a{Xs)ds - z 



biXs)ds)f{Yi),N,^n) 



wherefrom we obtain the basic decomposition of regular and singular part. The 
probability of taking jumps up to time 1 and starting at i is — which leads 
to the expression for the singular part of ii^'^ in 1)2. 9|) by the definition of /. Again 
by the definition of / , formula (|2.4|) and formula (|2.8|l we obtain that, for n > 1, 

£;(exp( f a{X,)ds - z f b{Xs)ds)f{Yi),Ni - n) 
Jo Jo 

= ^£;(exp(/" a{Xs)ds~z( b{X,)ds)f{Yi)\pn = l)P[Pn=l) 

= YEp^{cmI a{Xs\p-i^^^)ds-z f 6(X.|p-i(^))ds)/(Yi|^-i(^)))P(p„=.7) 

= V / exp(-z2;)— 77''^i'--''^"'''^i'''^i'-'''^"'''^i(da;)ord(7) 
^-^ In n'. 

7ec'„ 

holds true. We have applied that / equals de*-''^^-' ord(7) on the set ^^^(7), since 
this set consists of loops with favorable path 7 and Zi takes precisely the value 
ord(7). In addition de^'^-^^ P{p-\j)) = 
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We have to show that the measure ri°'ii'----'^i7i-''^ii''>ti,---,b-,n-b-,i (^d^) has a density, 
which can be seen from the fact that the intersection of A„ and the set {X]r=i + 
byitn+i = a;} depends in a differentiable way on x. The sum starts with n = 2 since 
there are no loops with only one jump. □ 

Example 1. We illustrate the stochastic approach by the case d = 2. Since a loop 
7 S C„ only appears if n is even and has the form 121 ... or 212 . . . , the con- 
tributions in the above series are necessarily non-negative: indeed for a hermitian 
2x2 matrix A we must have that ord(7) > for all loops 7 and the measures rj 
are non-negative either. Hence the density is non-negative. Again we note that the 
validity of the BMV- conjecture for d = 2 is well-known (see e.g. j2jj. 

Remark 5. Theorem\^can also be derived from the well-known Dyson expansion 
(see for instance j3]/ Still we believe that the stochastic reasoning underlying our 
proof has special merits: it leads us to a probabilistic and combinatorial point of 
view (see Section 3 "Stochastic Semantics"). 

Remark 6. We have formulated Theorem\^for Hermitian matrices A as this is 
presently our natural framework. But it is clear that it may as well be formulated 
for general d x d matrices A. 

3. Stochastic Semantics 

From now on we assume d — 3 and 61 = 0, we shall write Oi = an for i — 1,2, 3. 
In particular all matrices will be real from now on. From the point of view of 
stochastic processes we now have a dynamic picture of the problem to calculate 
the measure fj,^'^: we consider paths with values in the set {1,2,3}, which are 
favorable in the sense that two neighboring elements are different and the last 
element is different from the first one. The combinatorics of these paths leads 
us to a particular way to sum up the series (|2.1()|l . We think about trajectories 
dynamically as loops on the vertices {ei, 62, 63} of the 2-simplex. The total times 
fi, which the trajectory stays in state i during time [0, 1], form an element of the 
2-simplex. Only if X]i=i bi^i G [x,x + dx], the trajectory contributes to the density 
i;^'^{x)dx. 

We now fix < 63 < 62 and x e]0, 63 [. Due to the following choice of parameters 
we choose the unsual convention 63 < 62. The intersection of A2 with X]i=i bi£,i = x 
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will be parametrized by 

i ((1 - X2) + t{x2 - X3),X2{1 - t),X3t) 

with real numbers < X2 < X3 < 1 and t £ [0, 1]. We shall denote this line segment 
by ^^2,2:3 g^j^(j obtain the relations 

522:2 = X, 
bsxs = X. 

In particular we observe that - for given x - the numbers 62, 63 and X2, X3 determine 
each other. In order to obtain X2 < x^ we have been choosing 63 < 62- 

We apply the notions of the previous section. The characteristic char(7) = 
(^1(7), fc2(7), fc3(7)) of a path 7 e C is the number of visits in the points 1,2,3. 
Clearly fci + fc2 + fc3 — n. We shall observe in the following Lemma that in dimension 
3 the characteristic already determines the number of jumps between 1 — 2,1 — 3 
and 2 — 3, denoted by /12, ^13 and /23- These quantities are defined via 

'12(7) '13(7) '23(7) — n A( \ 

a^2 '^13 Ct23 '^7172*^7273 ■■■ '^7n-l7n'^7n7l — ora^7j 

and the numbers hjil) of jumps between i and j only depend on char (^7) = 
(A:i(7),fc2(7),fc3(7)) for7e C„. 

Lemma 1. The characteristic char(7) = (fci(7), fc2(7), ^13(7)) of a path 7 G C and 
the powers (^12(7), ^13(7), ^23(7)) of the order oid^'j) are in one-to-one relation. By 
abuse of notation we may therefore write hjij) — Zij(fci(7), fc2(7), ^3(7))- 

Proof. We take formula H2.7I) and solve it for lij given char(7), we obtain 

^12 = fci + fc2 - fcs 

^13 = fci + fcs - fc2 
^23 = ^2 + fcs - ki, 

which yields the result. □ 

Next we calculate in our particular setting (recall that d — 3 and bi — 0) 
explicitly the density of 7^a^^,...,a^„,a^j;b^^,...,&T,„,b^^ 
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Lemma 2. For ki, k2, fcs > define a probability density f on A2 (this time with 
respect to uniform distribution of total mass ^) given through 

(3.1) /(a, 6, 6) = /3(fci, fc2, fc3)^i'^"'C2'"'e3^"' exp(aiCi + 026 + ^36), 
where 

(fcl + ^2 + fcg - 1)! 



P{ki,k2, fcs) 



(fcl - 1)!(A;2 - l)!(fc3 - 1)! 
for ki > 1. We fix a path 7 G C with characteristic char(7) = (fci,fc2,fc3), n := 

fcl + fc2 + fc3, and define 7' = (71, ... , 7„, 71) g {1, . . . , and 

pr^ : A„ ^ A2 

through ^(7) := (pr^(ii, . . . , i„-(-i))i = Sj=i f'^^ * ~ lj2,3. Notice that the 

-l'j=i 

law of the real-valued random variable lu i— > b2^2{Pn{'-^)) + &3C3(Pn(^)) under is 
j^aT,j,...,a^„,a^j;b^^,...,7„,7H-^^-j _ n\(l)[j^x)dx (scc formula \2. 1 1\) Then the following 
assertions hold true: 

(1) Assume ki > 1, for i = 1, 2, 3, then the law of the random variable pr^ has 
density 

(3.2) (fcl + fc2 + k,)^fi^,,^2, 6) = rjf^/(a, 6, 6) 

with respect to the measure ^X2 on A2. Notice that the state appearing in 
71 is counted twice in the density. 

(2) For K>l,i = l, 2, 3, and x e ]0, b^l 



(3.3) 0(7, x) 




(n - 1)! V ^2^3 Jo 
(3.4) /(I - 2:2) + t(a:;2 - X3),X2{1 ~ t),xs,t)^x2Xzdt 

at < X < b^, where the cases in {} pertain to 71 — 1,2,3. Notice the 
relations b2X2 = ^32:3 = x. 

(3) Assume fci = 0, and x e]0, 63[, then 0(7,0;) = for all n>2. 

(4) Assume fc2 — 0, and x e]0, b^l, then 
(3.5) 



0(7, x) 



1 (fci + fc3-l)! I exp(ai(l ~ 3:3) + Q3X3) "tl"''"' »/7i = 1 



(n-l)!(fci-l)!(fc3-l)! I exp(ai(l-.3) + a3a;3)^i^^Cf^*/7i=3 
and n is necessarily even. 



fclf)3 
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(5) Assume ks = 0, and x s]0, 63 [, then 
(3.6) 



1 (A:i + fc2-l)! J exp(ai (1 - x^) + 02X2) ^'""ll"''" »/ 7i = 1 



(n - 1)! (fci - l)!(fc2 - 1)! ] exp(ai(l - .2) + a2.2)^i^^S-^ ^/7i = 2 



anrf n is necessarily even. 



Proof. Fix 7 e C„ and let 7' = (71, . . . ,7„,7i). We first set aj = for i = 
1,2,3. By direct computation we verify that now / indeed defines a probability 
measure on A2, hence the norming factor is correct (the actual form of / stems from 
pushing forward with pr^ and simply observing that a sum of independent uniformly 
distributed variables leads to a /J-distribution). In the chart 7ri2 (projection from 
A2 on the first two components in M^) the volume element ^X2{d^) equals d^id^2 

on{(a,6)eM2, a,6>0, el+6<l}• 



J/ /(a,6,6)A2(do = /3(fci,fc2,fc3) / / '^i'^-'d'-'(i-a-6)'='-'rfa'^6 

^ JA2 Jo Jo 

- Bik. ko ,j (fc^-l)Kfci + /c3-l)!(fei-l)!(/c3-l)! _ ^ 

We continue now with general ttj. Calculating the formula of the density ^(7, x) at 
X e]0, 63 [ amounts to calculating the mass of pr^ passed by the line L^^.xa through 
variations in x. Fixing 62,63 we thus fix 0:2, a;3. The area of the quadrangle with 
corners at ei + Xi{ei — ei), ei + {xi + dxi){ei — e{), for i = 2, 3, with respect to the 
measure \\2{d£) - under a small variation rfa; of a; - is given by 

\-{xzdx2 + X2dxz) = —^xdx 
2 6362 

= T-—\/b2hX2X3dx 

03O2 



-r—rVx2X^dx. 
02 03 

Shrinking the side convjei +2:2(62 — ei), ei +2:3(63 — ei)} to an infitcsimal clement 
at the point ((1 — X2) + t{x2 — .2:3), 2:2(1 — t),xzt), for t G [0, 1], on L^2'^3 leads to 
the appropriate area element 
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Hence we can determine (j)n(^,x) through equation H2.11|l and formula H3.2|l evalu- 
ated at ((1 — X2) + t{x2 — xz),X2{l — t), xst), for t e [0, 1], 

^'7(^26 °Pn + ^36 e[x,x + dx]) 




-r-^^/x^dx- — ^—-^f3{ki,k2, fca) 
O263 

/(I - X2+ t{x2 - X3),X2{1 - t),X3t)dt. 

For the degenerate cases we perform the same program. We first calculate the 
density of the law of pr-y if one of the ki is zero, which is a density supported by one 
edge of the simplex A2. Assume A:3 = 0. With respect to the uniform distribution 
with total mass 1 on the edge convjei, 62} of A2 we obtain for ki,k2 > 1 



and similar for the other case. A small variation dx in x leads via = dxi for 
i = 2, 3 to the desired results. □ 

In order to write the above densities in a more compact way we shall apply the 
well-known formula 

[\it)e-' ^ 9{o) 
r(a) Jo 

as a I for any continuous function g : [0, 1] — > R. Hence we can apply (fe — 1)! = 
T{k) for fc > and obtain the following proposition: 

Lemma 3. For 7 G C and x g]0, 63[, we obtain in the sense of Gamma- functions 
0(7, x) 




{n - 1)! V 6263 Jo 

/(I - X2 +t{x2 - X3),X2{1 - t),X3t)y/X2X3dt, 

for char(7) — (fci,/c2,fc3), ki + k2 + ks — n, and ki > 0, where the cases in {} 
pertain to 71 = 1, 2, 3. 

Proof. For fcj > 1 there is nothing to prove. Assume now that we take the limit 
^2 i 0, hence 72 = 1 or 3, since the vertex 2 cannot be starting point. We introduce 
furthermore 

A := ai{x2 - X3) - a2X2 + 033:3 
fi := ai(l - X2) + 02X2 
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as in Remark d Hence the limit yields 



(fci + fc3 - 1)! r^J_ / (^'"^)+*("^'"^) if 71 = 1 



S (n-l)!(fci- 1)1(^3 -1)! V 6263 r(a) Jo \ ^ if 7i - 3 

((1 - X2) + t{x2 - X3) f'-\X2{1 - t)r-\x3tp-^^^dt 

exp(A + /i) (ki + h-iy. 1 /iiij / ^i^^^^J^ if 71 = 1 



(n-1)! (fci-l)!(A:3-l)!a;2V ^2^3 ] iiz^ll!!!^ if = 3 

^ exp(A + M) (fci + fc3-l)! / ^'~''Cbf'' if 71 = 1 
(n-1)! (fci-l)!(fc3-l)! ] (lzf3)!^if 3 ' 

since X2&2 = a^3fe3 = a;. Similarly for the third case. □ 

For the calculation of the BMV-measure ^'^'^ we can make an essential further 
simplification: it turns out that if we average over all paths 7 with fixed characteris- 
tic char(7) (and varying the first entry 71) formulas H3.3|l - H3.6() appear in a simpler 
form, which only depends on the characteristic. We define the density 

X{ki,k2,k3,x) := ——— , \ . j-, — , , V (t>{l,x), 

#{7 e C : char(7) = (fci, fe, ^3 } ^ 

char(7) = (fei ,k2 .fca) 

i.e. the average of the densities 0(7, x) where 7 ranges through the paths with fixed 
characteristic {ki,k2,k3). 

Lemma 4. We fix a path 7 G C with characteristic char(7) — (fci, fc2, fc3) for n > 2. 
Then the following assertion holds, 

(3.7) X{h,k2,k3,x) = , I f{l-X2+t{x2-X3),X2{l-t),X3t)y/x^dt 







62^3 " ! 

in the sense of Gamma-functions 



Remark 7. The exponential term in simplifies to e^^^^ with X = ai{x2 — 

X3) — 0,2X2 + 030:3 and /i = ai(l — X2) + 0,2X2- Hence we obtain, for a path with 
characteristic char(7) = (fci, /c2, fe), ki > 1, by the binomial theorem and the Beta 
integral that 



Xiki,k2,k3,x) 



(1 -xa )'-■!- 



nx 



EE 



— ^ A-^ /fci - 1\ /a;2 - a;3 V (fc3)L+r 



^ Ll\ r \l-x2 (fci-l)!(fc2 + fc3 + i + r-l)!' 

L>0 r— \ / \ / \ / \ / 
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holds true, or — hy using t — 1 — s — an alternate representation, 

X{ki,k2,h,x) 
fei-] 

EE 



L>0 r=0 



nx 



Al - A n+X (-A)^ ( X:i-X2 \ ' (fc2)L+r 

n \\~x^) (fci-l)!(fc2 + fc3 + i + r-l)!" 



r(r+fc) 



_ffere we apply the notion {k)r '■= r(fc) 

Proof of Lemma 3. For the proof we apply the representations of the densities 
(|3.fi(l . and the fact that among all paths 7 G C„ with characteristic char(7) — 
{ki, ^2, fca) the path with 71 = i appear with relative frequency hence absolutely 

#{7eC: char(7) = (fci,fc2,fc3)}- 

n 

times. We calculate the density x{ki, k2, ks, x) at x e]0, 63[. This leads for h > 1 
to 

Xiki,k2,k3,x) = ————^ V 0(7,3;) 

#{7 e C : char ^ 

(7) = {ki,k2, fcs)} char(7) = (fci,fe2,fc3) 

1 n~ f\ki{l-X2)+tix2-X3) ^ 



[n - 1)! V fe2&3 Jo ki 

^ k2 3:2(1 - _^ fc3 ^3^\ 

n /c2 n k^ J 

/(I - 2:2 + t(x2 - X3),X2(1 - t), X3t)^X2X3dt 

1 r' 



/(I - 2:2 + t(a;2 - a;3),a;2(l - t),X3t)^x2X3dt. 

n\ V 0263 Jo 

For fci = we conclude directly. For fc2 = we use 

„ , , 1 k, {k,+k3-l)\xl'-\l-X3f^ , , , 

x(fci,0,A:3,x) = — — ^ exp(ai 1 - 3:3) + 03x3 + 

(n— Ij! n A;i!(fc3 — Ij! 63 

, 1 + 



(n-1)! n (fci-l)!A;3! 63 

■ exp(ai(l - 0:3) + 030:3) 



1 (fci + fc3-l)! x)^-\\-X3f'-^ 



n\ (fci - l)!(fc3 - 1)! 63 
and analogously for ^3 = 0. □ 

For the case x g]63, 62 [ we shall apply the following parametrization 



t^{{\~ t)yi, (1 - yi) + t{yi - ys), tys) 
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for < j/i < 2/3 < 1 satisfying the relations 

hyi = b2-x 
ih - bz)yz = b2- X. 

This leads as in the proof of Lemma [SI to the volume element 

1 



y/b2{b2 - 63) 

under variations of x, hence the respective densities x satisfy the following relations: 
we fix a path 7 e C„ with characteristic char(7) = (fci, k2, fcs), ki + k2 + — n for 
n > 2, hence 
(3.8) 

X{ki,k2,k3,x) = J j—T T—^ / f{{l~t)yi,l-yi+t{yi~y3),ty3)^/yry3dt 

y 02(02 ~ 03) n\ Jq 

in the sense of Gamma- functions. 

Remark 8. Notice that the case x g]63, fo2[ deduced from the case x e]0, 62 [ by the 

permutation 1 < > 2 is performed. One replaces then X2 by yi, x^ by y^, performs 

the permutation for aij , and replaces 62 by b2and 63 by &2 ~ ^3- All the necessary 
relations maintain and the first case in full generality then implies the second one. 

4. Combinatorial Sums 

Our next goal is to represent ^^'^[x) := ipix) in the following way. By Remark 
IHlit suffices to consider the interval ]0, 63 [. 

Proposition 3. Suppose that 62 > ^'3. Then, for x e]0,63[, we have 

^(^) = X! X{ki,k2,k3,x) ord(7) 
-tec 

fc>l m>0 />0,i=m mod 2 L>0 

fc-l /, ,x / \r (2k+rn--l\ 

'k-l\ Xj,-X2\ \ 2 ) 



E 



Ir+L 



r J \ 1 — X3 J k\(k + m + r + L — ly. 



r=0 



J \ fc - 1 

0<.j<k.j—m mod 2 

X (ai2Vii)^''~' (ai3\/^)' (0-23^/^2X3) 



2 
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The proof of Proposition|3|is just a direct combination of Remark|Z| the following 
Lemmainiand the representations Zi2 = 2k~l, k2 = (2k—l+m)/2, and fca = {l+m)/2 
when ki = k, hs — I, and I23 — m are given. 

However, the representation of ^p{x) in Proposition |3| has to be transformed in 
a proper way to observe that it is non-negative. For this purpose we will further 
introduce the hypergeometric function F{a, b; c; z) and use certain hypergeometric 
identities in order to simplify the above representation. 

4.1. Counting paths on the triangle. 



Lemma 5. The number of paths ^ in C with ^1(7) = k, ^13(7) = /, ^23(7) = Tn 
and I = m mod 2 is given through 

2k + m ^ fk\l"-^ + k-l\fk-f 



(4.1) 



E 



fc- 1 



(}<j<.k.j=m mod 2 

If I ^ m mod 2, then the number of paths vanishes. 



1-] 

2 



V. 



Proof. From ^ we get that the generating function of ord(7) of all paths 7 with 
71 = 1 is given by 



7.71=1 







1 

023 


023 
1 




1 

-ai2 
-ai3 


-ai2 

1 

-023 


-O13 

-023 
1 



1 - 023 



1 — 2ai20i3a23 — 0^2 — O13 — O23 



Hence, if Pi{k, Z, to) denotes the number of paths 7 in C with ki{'-/) = k, /i3(7) = Z, 
^23(7) — m, and 71 = 1 we have 



^ Pi{k,l,m)x''a' 

k,l,7n 



I rn 



1 - 0^3 



" 1 - - a;(2ai3a23 + 0^3 + 1) ' 

From that we immediately get (if / = m mod 2) 

0<j<k,j=m mod 2^"'^^ /\2/ 

Finally, if we denote by P(fc, I, m) the total number of paths 7 in C with ki (7) = fc, 
^13(7) — I, and ^23(7) = ™ then 

—Pi{k,l,m) = —P{k,l,m), 
k n 



where rt = 2fc -|- m = fci + fc2 + fc2 • This proves H4.ll) . 



□ 
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4.2. Hypergeometric identities. The hypergeometric function F{a, b; c; z) is de- 
fined (for complex \z\ < 1) by 

where (x)„ = r(a; + n)/T(x) = x(x + I) ■ ■ ■ (x + n — 1) denote the rising factorials. 
There are lots of identities (see Chapter 15]) for these kinds of functions. Some of 
them will be used in the sequel. For example one has Euler's integral representation 

F{a, b; c; z) = Jf'^ f\l - zt)-'^t''-\l - t)^-'-^ dt 
T{b)T{c-b) Jo 

if \z\ < 1 and 5R(c) > 5R(6) > 0. Furthermore, it was already known to Gauss that 

, . ^, , ^ r(c)r(c — a — 6) 

4.2 F{a,b;c;l)^ ) ' ^ -( 

i (c — a)l (c — b) 

if 5R(c-a-6) > 0. 

We start with a lemma, where we use the identity 

z 



(4.3) F(a,5;c;z) = (1 - |^a,c- 6;c, ^ ^ 

Lemma 6. Suppose that j = m mod 2. Then 

2^ E fc^f^^^) c'i.-- 



1^3:1—3 mod 2 



(fc-i)! (C 

V 



w/iere = 2CD/{C^ + D^). 

Proof. We note that the left hand side of the above equation can be represented as 

' 2k + m — 



xF[-ik-j),- 



2 

2k + m- j \ f 2k + m- j \ C 



-i2 



2 7 ' V 2 ' ■ 7 ' D2 

and the right hand side as 

' 2k + m — j 



{C^ + D^f-^{2CDy 



2 



X F -(fc - j), -r; - + r - 1 ^ ■ 



C2 + D2 



By using H4.3|) with 

f2k + m-j \ (2k + m-j 
a = ~{k-j), b^- I 11, c = - ( + r- 1 

and z = —C^/D^ we directly get a proof of the lemma. □ 
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4.3. Further Hypergeometric Identities. In this section we present a proof 
of rather strange identities that seem to be new in the context of hypergeometric 
series. 
We set 

Mk; V, := E E { ) t ^^"^ - 

"Vn'^VJ/ (m + l)fc+2r-l m! 

where r a is non- negative integer. 

Lemma 7. We have 

Ao{k;v,i) = {l + vfe^ 

g wiw^i' i}^T (^)' 

and 



i^^m + r-\)\{k-2l)V 

-io (r-1)! \ 2 J 2' 2 (1 + ' 

where r a is positive integer. 

Remark 9. Note that the right hand sides of these identities are non-negative if 
\v\ < 1. Hence, we have Ar{k;v,^) > 0. 

In fact, we are more interested, in sums of the form 

Ar{k;v,0 = ^{Ar{k;v,0 + Ar{k;-v,-0) 



k\ ^ [ 2 )k+r-l 



^0 0<J<k, l^m mod 2 V^V + l).+2r-l m! 

Since Ar{k; v,^) >0 and Ar{k; —v, — ^) > O(for |v| < 1^ we also have Ar{k; v,^) >0 
and the representations 

fk\l-V^ ,k-2r.f k-2 k - 3 ^ {1 - V^){1 - S^)\ ^ 
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Ar{k-V,i,) 



_i (r-1)! V 2 y V 2' 2 ' ' (1 + 5^)2 

where r a is positive integer. 

Proof. ^ We prove first the case of positive r. Botfi sides of the identity are power 
series in v and ^. Thus, it is sufficient to compare coefficients. The coefficient of 
v^£,™/ml of the right hand side is given by 



^^^e\ii + r-iy.{k~2iy. 



X 

i>Q 



By applying the substitution s = \/t, integrating the corresponding Beta integrals 
and rewriting the sum over £ in hypergeometric notation we thus get 
fc! 



^J^ie + r-mk-2£y. 



l>0 

X 



i>0 

fc! 



^^^£!(£ + r-l)!(fc-2£)! 



X 

i>0 



A 1 /fc-2Ar(m/2 + j/2-i + l/2)r(^ + r) 



^ij 22*+'- \ j -2iJ T{£ + r + m/2 + j/2 -i + 1/2) 
{-ly{l+i)k-^ 



E 22.+. _ 2z)!(fc - + z + § + ^y+r 



"^U"2'2 + 2-2'2 + 2 + Y+^'l 
Next we use formula (|4.2() and obtain (after rewriting the remaining sum in hyper- 
geometric notation) 

fc^^(i + i + f)^(-i + fc + f + r)^^ J 1 J 1 J 



j7 2'T(| + f +r)r(i + | + f +r) V 2' 2 2' 2 2 2' 



^This nice proof was pointed out to us by Christian Krattenthaler and is considerably easier 
than our first one. 
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In order to avoid difficulties with zero-cancellations we interprete this sum as a 
limit, use again formula H4.2|l and obtain (after some algebra) 

lir, A\ r(^ + i + f)r(-§ + fc+f +r) / J l^l.l^l^rn 
e-oVi^Td + f +r)r(i + | + f +r) V 2'2 2' 2 2 2 ^^' 

k\ 2'^-+^+2^-2p (_i + _^ J. _^ m) p _ 2£ + 1) sin(7r(2e - m)) 



lim , , 

\i) (fc + m + 2r)!r (l - 1 + f - e) sin (vr (i ~ § - f + e)) sin (tt (§ - f + e)) 
Now note that the limit of the sin-terms is always 2. Hence, we finally obtain 

ofe+r-l ( m-j+2 \ 
V 2 )k+r-l 

(m + l)fc+2r-l 

as proposed. 

The proof for the case r = runs along similar lines. The only difference is the 
singular term (^) in front. However, after integrating the Beta integrals we can 
rewrite the corresponding sum as 
k\ ^ fc! 



X 

i>0 

fc! 



1 /fc-2^-2\ r(TO/2 + j/2-z-|-l/2)r(^-f 1) 



'\k y 22^+2 V j-2i ) r(^ + m/2-hj72-i-h3/2) 

E 



^^_-^^!(^+l)!(fc-2^-2)! 



X 

i>0 



1 /fc-2^-2\ r(TO/2 + j72-i + l/2)r(^+l) 



k j2^'^+^\ j-2i J r{£ + m/2 + j/2-i + 3/2) 
(-l)'(l +i - 2i)fc_j+2, fj k 1 j kljm^ 



|^22«(fc-j)m(i-z + l + f), ^2 2'2 2 2'2 2 2' 
and proceed as above. □ 

Lemma 8. Set 

(4.4) 

T(fc, p; ^) := E g g (^ + l)Vr"" W^^- ^ 

(4.5) S{k, r, p; C) := ^(-l)'^'"-^ " ^) T{k, r,p + T; v, 0- 
Then 

(4.6) Tik,r,p;v,0^J2[ ^ ) S{k,r, p + t;v,0 



24 MICHAEL DRMOTA, WALTER SCHACHERMAYER, JOSEF TEICHMANN 

and 

(4.7) Sik, r, p; 6 = E "/) ^^A.,,ik - 0- 

In particular we have S{k, r, p; w, ^) > and T(fc, r, p; w, ^) > if \v\ < 1. 

Remark 10. If we set T{k,r, p;v,£,) = ''i + "^(^j '^i P! ~0) '^"'^ 

S{k, r, p; V, ^) = ^{S{k, r, p; v, £,) + S{k, r, p; —v, — C)) then we have (of course) cor- 
responding representations in terms of Ar{k;v,^) and also S{k,r, p;v,^) > and 
fik,r,p;v,O>0 tf\v\<l. 

Proof. First note that H4.5|l and H4.6|l are equivalent. Thus, it remains to prove 
(|4.7() or equivalently 

T{k,r, p;v,Cj 

"^fr-p\ kl ^ /r-p-r\ (2a)! , 

By expanding both sides with respect to (J^ /m\ this identitiy is equivalent to 

k\ ^ V 2 )k+r-p-l [k - jy. 



jj {m+l)k+r-i {k-j-py. 

(r-p\ kl (r-p- r\ (2a)! (k-p-r\ 2^+-^ 
' ^\ ^ ) {k- p-t)\\ 2a ) 2-a\\ j ) (m + l)fc+2a+p+r-i 



T,a>0 



By rewriting the sum over r of the right hand side in hypergeometric notation and 
by using ()4.2|l we get 



a>o •^''^^ ~ P~ -^^'(^ ^ ^ ^ 2a)!a!(m + l)k+2a+p-i 

X (-(r - /? - 2a), -(fc - /? - j); TO + fc + 2a + p; 1) 



E ii 



a>o -^'^^ ^ ^ ^ ^ ^ ^ 2a)!a!(TO + l)fe+2a+p-i 
r(m + fc + 2a + p)r(TO + 2fc + r - p - j) 
r(TO + fc + r)T{m. + 2fc + 2a - j) 

Next this sum can be also written in hypergeometric notation. Further, a second 
use of 1)4. 2|l and some simplifications (using the duplication formula of the Gamma 



A HYPER-GEOMETRIC APPROACH TO THE BMV-CONJECTURE 



25 



functions) yield 

^I2fe-i (^inzJ±l') m\T{m + 2k + r - p - j) 
j\{k - p- j)!r(m + k + r)r(m + 2fc - j) 

„ ( r — p r — p — I rn + 2fc — 7 + 1 

X F ; : 1 

V 2 ' 2 ' 2 

k\ {k - j)!2'=-i r(m + 2fc + r - p - j) 



{k- p- i)\{m + l)fe+^„ir(m + 2fc - j) 

r(fc + ^-| + i)r(fc + ^-| + r-p-i) 
r(fc + ^ — | + | — | — i)r(fc + ^ — | + | — | — 1) 

jj (m + l)fc+r-i (fc-j-p)! 

as proposed. □ 



5. Proof of Theorem ^ 

First we use the results of the previous section to obtain another representation 
for il>{x). 



Lemma 9. Suppose that 62 > ^3 and set 

A12 = ai2V^ = «12 
^13 = ai3V^ = «13 

2A12A13 

^12 + ^13 



X 
X 



£. = a23\/X2X3, 
X3 - X2 

W2 



2{l~X3) 



A12 + A\^ 

Then for x g]0, 63 [ we have 

^ ■^^ k>l ^ ' ■r=0 ^ / L>0 



r+L 

X 

p=0 



E f'^^^)5(fc,r + L,p;«,e)c.''. 
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Proof. With help of Proposition 13 and Lemma El we get 



x(l - X3) k\ ^ \ r / L! 

^ -^^ fc>l r-O ^ ^ L>0 



j7V fc-i 



m>0 0<j<k,j=m mod 2 
xVr (^"^ + gfc+r+L-p-l / 2fc + m - j 



Efc!(fcllvE( ^ ) ^2E 



^ k>l ^ ' r^O ^ ' L>0 

'"+-^ / I r\ / /I \ P 



■ s ofc+r+L-p-1 / m-j+2 \ 



X 

p=0 



j7 (m + l)fe+,.+L-i (fc-j-p)! m! 



xE E 

m>0 0<.j<k.j—m mod 2 

i-^ /fc _ A ^ (-A)^ 
:r(l-X3)^fc!(fc-l)! ;^^l r j ^ L! 



L>0 



X 

p=0 



Finally by using H4.6|l we directly derive the proposed representation. □ 

Note that |t;| < 1. Thus this lemma shows that ip{x) > if > and A < or 
equivalently \Ai2\ > l^ial and 01(62 — 63) + 02^3 — a3&2 > 0. This is satisfied by 
assumptions of Theorem^ Hence we have proved Theorem^for x g]0, 63[. The 
case X (z]b3,b2[ is done by exchanging indices 1 and 2 and 63 by 62 — ^3 (compare 
with Remark 121). 

6. Appendix 1: the Feynman-Kac formula 

We shall work with a continuous-time Markov process on the state space S — 
C X {1, . . . , d} associated to the off-diagonal elements of a d x d matrix A G Md{C) 
with complex entries. The non-zero diagonal elements of A and the "potential" B 
will appear in exponential functionals of the process. We denote C-valued functions 
on the state space 5* by /. Given a matrix A g Afd(C) with zero diagonal, we 
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associate a generator A of a pure-jump type process on S, namely 



Af{C,i) = J2{f{Caij,j)-f{C,i)). 



The resulting 5-valued Markov process is denoted by Y^^'^^ := {Z^^'^\x^^'^^) with 
initial value {(, i) G S and has the following properties. 

Lemma 10. The projection of (Y}''''^)t>o =: {zi''''\ X^''''^) t>o onto the second 
component will be denoted by {XI), since it does not depend on C. G C, and equals 
in distribution the {1, . . . , d} -valued Markov process associated to the matrix G with 
entries gtj = 1 for j ^ i and gu = 1 — d for i,j = l,...,d. 

Proof. Take a function / on S, which does not depend on the first component C, 
then we have 

^/W = E(/(i)-/W) 

= {Gf)u 

where we can identify / with a vector in C''. □ 
Lemma 11. Let f^:Cx {1, . . . , rf} — > C be the function 

for {C,i) G S and r G (which may be identified with a linear map from C"^ to 
C{S,C)), then 

Af = f(A-(d-l)U)r_ 

Proof. Direct calculation of the generator A on functions /'". □ 

Furthermore given a real valued function on S of the form 

V''{<:,i) = bi 

for (C,*) G 5* (playing the role of a "potential") and b G C'', then we can define the 
multiplication operator (denoted again by V^) and we obtain 

■^6 j^r f^^ 

where B denotes the diagonal matrix with entries b\,...,bd € C. 
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The main assertion of this section is the Feynman-Kac formula for the family 
of Markov process (Yf '^'^^)t>o for (C,i) G S. Since we do not want to go into the 
theory of Feller semigroups and maximal domains, we formulate the Feynman-Kac 
Theorem in a special case, i.e. on a finite dimensional domain of definition. 



Theorem 3. Let A be any d x d matrix with entries in C and diagonal elements 
fli, . . . , Gd, let be the multiplication operator for b E C^' as defined above, then 
for z €E C the function, 

= E (^exp(£ a{Xl)ds - z j\{Xl)ds)f^ {Y,^^''^)^ 

solves the following differential equation on C{S,C) for initial value f^, r G , 
and all < > 0, 



(6.1) ^ut = Aut + Vut - zV^ut 

at 

(6.2) = /• 



On the other hand, for initial value f = f^, the solution to the differential equation 
is given by 

t ^ e^tid-l) px.p{t[A-zB))r ^ 



henc 



E(^exp{J^ a{Xl)ds-z b{Xl)ds)r{Y}^''^)^ ^ e-*('^-i)(exp(i(A - zB))r), 



for all z G 



Proof. The proof is done by the Markov property for the process (y/''''-')f>o with 
(C, i) G S: first we show that the following semigroup property holds: 

ul^+.^iCi) = E{eM r a{Xi)ds - z r b{Xl)ds)ul{Y,f'^)) 
Jo Jo 
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for ti,t2 > 0. Indeed, the right hand side can be written by the Markov property 

i?(exp( a{X:)ds - z b{X:)dsWMf% 
Jo Jo 

/ rt2 rt2 pti 

= E exp( / a{Xl)ds ~z j b{Xl)ds)E{exp{ / a{Xl)ds- 
\ Jo Jo Ja 

('exp( / a{Xl)ds -z [ b{Xl)ds)E{e^p{ [ a{Xl)ds- 
\ Jq Jo Jti 



E 



tl+t2 



-z I b{Xl)ds)r{Yl^fj\T,J 
Jtl 

/•tl+*2 /•«l+t2 



= i?(exp( / a{Xl)ds - z I b{Xl)ds)) f {Y,%1)) 
Jo Jo 

~ fti+t2 (Cj 

for aU ti,t2 > 0, (C, «) e S* and r G C^. Furthermore r i-^ u^C, i) is obviously hnear 
in r for fixed (C, i) ^ S and t > 0. Hence it is sufiicient to calculate the derivative 
at t = 0. This leads to 

^\t=oul{C,i) = air,(^ - zbir,C + Au^iCi) 

Hence (uj )t>o.reC'' defines a strongly continuous semifiow on the space {f^,r g 
C*} with generator A + — zV^, which a fortiori coincides with the semifiow 

t ^ e"*(''"l) yoxp(t(A-2S))r^ □ 

7. Appendix 2: the first conjecture 

We provide a counter-example for the original version of the BMV-conjecture as 
stated in P]. We work with the notations of Section 3. Given the matrices 

A = 

and 

B = 

we define the signed measure ^f'^{dx) as inverse Laplace transform of z i— s- (ei, exp(^ — zB)ei) 
In the notations of Section 3 we have 6i = 63 = and 62 = 1. We calculate the 
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A,B 



sign of the absolutely continuous part of ' asymptotically in e and show that we 
obtain a negative sign for small e. We apply the formulas in the sense of Remark 
IHl We notice that = j/3 = 1 — a; for < a: < 1, which leads to the formula 



(j){ki,k2,ki,x) 



/((l-i)(l-x),x,t(l-x)) 



("-l)Uo 

for trajectories 7 and 71 — 1, hence by Corollary^ 



E 

■yec 

n>2,7i = l 



(l>{ki, fc2, fcs, x)a[^^a[^^a''^.^ 



We only have to calculate the following cases up to order e'^, where # denotes 
the number of possible paths 7 with given lij , where we apply Lemma El hence 
ki = h2±h3 > 1 and 2li2 + hs + ha < 4. This leads to the following table, 



^12 


ll3 


^23 


# 


2 








Pi(l,0,0)-1 


1 


1 


1 


Pi(l,l,l)-2 





2 


2 


Pi(l,2,2) = l 





4 





Pi(2,4,0) = l 



associated to the paths 121; 1321, 1231; 13231; 13131. The fourth path leads to 
vanishing (j) since the vertex 2 is not visited (see Lemma El 3. in the appropriate 
translation as in Remark jSJ. Hence we obtain up to order the following density 
for the absolutely continuous part 



(x) = (1 - x)C^r - 2(1 - x)h''^'- + (1 - xfe'l + Oie% 



consequently 



Ui^t'^'ix) 



3(1 -x)- 6(1 - x)2 + 2(1 - xf + ©(e) 



where we obtain a negative sign if x is near 0. 

We could also calculate in the original stochastic language by analysing jump 
densities and calculating with conditional expectations, which is instructive in the 
given example. We apply the formula of Corollary ^ and the definitions of Section 
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2, 

e-^{exp{A- zB)ei)i = E (^expij^ a{X,)ds-zj^ h{X s)ds) f {Yi)\Xa = I 
- E (exp(-z^ l{x,=2}ds)r{Y^)\Xo - 1^ 
Ep^{eM-z f l{x.^2}ds)\^i - l)P{Pn\l)\li = l)«i^i^«i^l4f + 0{ 

71 = 1 

2!l2+il3+i23<4 

where (^s)o<s<i is the standard Poisson process with intensity (d— 1) on the state 
space {1,2,3}. Now we can evaluate the sum directly: 

• the probability for a path configuration 7 with n jumps is , see formula 
12.81 conditioned on 71 = 1. The duration under the condition that the path 
configuration is 7 is uniformly distributed on the simplex A„ (see Section 
2 for the precise statement). Rescaling due to the factor e^^ on both sides 
leads to the precise formulas. 

• the loop 1 leads to the singular part of /^-^ ' . 

• the loop 121 yields the first non-trivial term in ' . This non-trivial con- 
tribution is the distribution of the total time in state 2 under the condition 
that the trajectory is 121. The density of the projection from A2, if we fix 
the uniform distribution on the 2-simplex, 

PiT2-Ti <x)^ ( d\2{tiMM) 

JA2,t2<X 
l-X <.l-t2 ^2 

= 2 / / dtidt2 = 2{x ). 

^0 Jo 2 

The rescaled probability of 121 is ^. The density of the second projection 

is 2(1 — x), hence the first contribution (1 — a;)a^2- 

• the loops 1321, 1231 yield the second non-trivial contribution, each of it 
coming from the density of one projection from the 3-simplex A3, which 
yields in the first case 



PiTs -T2<x)= [ d\^{h,t2, h, U) 



I A3,t2<X 
/•X rl—t2 /•l — t3—t2 

6 / / / dtidt3dt2 
Jo Jo Jo 
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and in the equal second case P{T2 — Ti < .t) = 1 — (1 — x)'^. The rescaled 
probabihty of each path 1321, 1231 is ^ , hence for the two paths we obtain 
the contribution (1 — x)^ai2ai3a23. 
• the loop 13231 yields similarly via the density of the second projection 
the function 4(1 — x)^ with rescaled probability the contribution (1 — 
x)^al^ai2. 
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